arXiv: 1502.03928v2 [cond-mat.supr-con] 1 Mar 2015 


Theory for superconductivity in alkali chromium arsenides A2Cr3As3 (A=K,Rb,Cs) 


Yi Zhou, 1,2 Chao Cao, 3 and Fu-Chun Zhang 1,2 

1 Department of Physics, Zhejiang University, Hangzhou 310027, China 
2 Collaborative Innovation Center of Advanced Microstructures, Nanjing 210093, China 
3 Department of Physics, Hanqzhou Normal University, Hanqzhou 310036, China 

(Dated: March 3, 2015) 

We propose an extended Hubbard model with three molecular orbitals on a hexagonal lattice 
with D^h symmetry to study recently discovered superconductivity in A2Cr3As 3 (A=K,Rb,Cs). 
Effective pairing interactions from paramagnon fluctuations are derived within the random phase 
approximation, and are found to be most attractive in spin triplet channels. At small Hubbard U 
and moderate Hund’s coupling, the pairing arises from 3-dimensional (3D) 7 band and has a spatial 
symmetry f y t 3 X 2 _ y 2 \, which gives line nodes in the gap function. At large U, a fully gapped p- wave 
state, p z z dominates at the quasi-ID a-band. 

PACS numbers: 74.20.-z; 74.20.Mn; 74.20.Rp; 74.70.-b 


Recently, CrAs based superconductors have attracted 
much attention. CrAs itself is a 3D antiferromag- 
net (AFM), which becomes superconducting (SC) un¬ 
der modest pressure with T c ~ 2K0. Subsequently, 
a new family of quasi-ID superconductors A2Cr 3 As3 
(A=K,Rb,Cs) has been discovered at ambient pressure 
with T c up to 6.1Kifl. The key building block of 
A 2 Cr 3 As3 is the ID [(CrsAss) 2- ]^ double-walled sub¬ 
nanotubes (Fig. [TJ, which are separated by columns 
of K + ions, in contrast to the layered iron-pnictide and 
copper-oxide high Tc superconductors. 

This new family exhibits interesting properties in both 
the normal and SC states. In the normal state, the resis¬ 
tivity follows p(T) = po + AT in a wide temperature 
region, different from the usual Fermi liquid behavior 
Po+AT 2 ^ H. NMR measurements on K2Cr 3 As 3 showed 
a power-law temperature dependence 1/Ti ~ T 0 ' 75 above 
T c , which is neither 1/Ti ~ T for a Fermi liquid nor 
Curie-Weiss behavior 1/T\T ~ C/(T + 9) for an AFMjj|. 
Below T c , the electronic contribution to the specific heat 
C e {T ) deviates from the BCS senario, and the extrapo¬ 
lated upper critical field H c 2 exceeds the Pauli limit [2j- 
Qj- The Hebel-Slichter coherence peak of 1/Ti is ab- 
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FIG. 1: Crystal structure of A.2Cr3As3. (a) Top view and (b) 
the (Cr6As6)oo sub-nanotube. In panel (a), larger and smaller 
atoms are at z=0.5 and z=0.0, respectively. In both panels, 
green balls indicate alkaline atoms, red balls arsenic atoms, 
and blue balls chromium atoms. 


sent in K^CrsAss NMR measurement [Bj]. London pen¬ 
etration depth measurement for K2Cr3As 3 shows linear 
temperature dependence, AA(T) ~ T, at temperatures 
T <T C 0. All these experiments are very difficult to ex¬ 
plain within electron-phonon coupling mechanism, and 
suggest unconventional nature of superconductivity. 

The electronic structure of K 2 Cr 3 As 3 has been inves¬ 
tigated by Jiang et al.Q using density functional theory 
(DFT), which is confirmed by later calculation Q. There 
are three energy bands at the Fermi level: two quasi- 
ID a- and /3-bands with flat Fermi surfaces, and a 3D 
7-band. It is natural to ask the question, which band 
is responsible for the superconductivity? On one hand, 
the linear temperature dependent resistivity and the 
power-law 1/Ti in NMR indicate a quasi-ID Tomonaga- 
Luttinger liquid; on the other hand, upper critical field 
measurement H c 2 implies a 3D superconductor In this 
Letter, we shall carry out theoretical study to address 
this issue. We shall focus on K 2 Cr 3 As3 and the theory 
may be extended to other alkali chromium arsenides. 

We begin with constructing the Hamiltonian by using 
symmetries. The space group for A 2 Cr 3 As3 lattices is 
P6 to 2 and the corresponding point group is D$h. 0-0 
We shall also assume that the time reversal symmetry 
remains unbroken, and consider a system described by 
the Hamiltonian 


H — Hq + Hi n t, 

where Hq is the non-interacting part and Hi nt is the in¬ 
teracting Hamiltonian. 

Tight-binding model. We assume Hq to be given by 
a tight-binding Hamiltonian. For K2Cr3As 3 , there are 
three Fermi surfaces corresponding to the a, /3 and 7 
bands. Minimally, we need three orbitals per unit cell 
(per CrgAs6 cluster). From the DFT calculation, there 
are three low energy molecular orbitals. Two of them 
belong to 2D irreducible representation E' of D^h group, 
and the other one is in ID representation A[. We de¬ 
note two E' states as |1) and |2) and the AI X state as |3). 
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Neglecting spin-orbit coupling, the tight-binding Hamil¬ 
tonian Hq can be constructed from these three orbitals, 

Ho = 'y ' C kms £kmnCkns) (1) 

k mns 

where m, n = 1,2,3 are the three molecular orbitals, 
s =t4 is the spin index, c kms (ckms) creates an m-orbital 
electron with spin s. In the basis {|1), |2), |3)}, (kmn can 
be written in terms of matrix form 

8 

£k — 5 k Ap (2) 

T = 0 

where Ai_§ are Gell-Mann matrices. 

We now use the symmetry to analyse £ k , r = 1, ■ ■ ■ , 8 . 
Ai_g transfer as irreducible representations under D^h 
symmetry operations. Namely, Ao and As transfer as 
A[, A 2 transfers as A! 2 , (Ai,A 3 ), (A 4 ,Ae) and (A 5 ,Ay) 
transfer as E' respectively. Hq should be invariant un¬ 
der all the Doh symmetry operations, thereby belongs 
to representation A\ . Therefore, 5£ can be determined 
using the Clebsch-Gordan coefficients. 0 Since the ab- 
plane lattice constant a = 9.98A is much larger than 
that along the c-axis c = 4.23A, we will only keep the 
hopping terms on a 6 -plane up to the first nearest neigh¬ 
bor (NN) bonds and those along c-axis up to the second 
NN bonds. To do this, we set the Bravais lattice ba¬ 
sis a = (§,— : ^, 0 ),b = = ( 0 , 0 ,c) (see 

Fig.l), and introduce k a = k ■ a, kb = k ■ b, k c = k • c, and 
some harmonic functions on the ab plane, 

s x 2 +y 2 (k) = cos k a + cos kb + cos(k a + kb), 
p x { k) = 2 sin(fc a + kb) + sin k a + sin kb, 
p v { k) = -\/3(sin kb — sin k a ), (3) 

d x 2 _ y 2 (k) = cosfc a + cosfcb — 2cos(fc a + kb), 
d xy { k) = — Vs(coskb — cosfc 0 ). 

It is easy to verify that s x 2 +y 2 oc 1 — j(fc^ + k y ), p x oc k x , 
p y oc k y , d x 2_ y 2 oc k x — k y and d xy oc k x k y at small ka. 
Thus, we have £ k = £ k = 5k = 0, and other nonzero £ k ’s, 

^k = 2 <i cosfc c + 2t4COs2fc c -(-2t 2 s a; 2 +y 2(k), 

5 k = 2 1\ cos k c + 2t\ cos 2 k c + 2 t' 2 s x 2 +y 2 (k), 

£k = -Ap y (k) - Bd xy (k), (4) 

5k = Ap x (k) - Bd x 2_ y 2{k), 

5k = Cp x (k) - Dd x 2_ y 2(k), 

5 k = Cp y (k) + Dd xy (k), 


where A = (£3 + to cos k c ),B = (£3 + £ 5 cosfc c ),C = 
(£3 + £5 cos kc),D = (£3 + £5 cos k c ). By fitting the DFT 
band structure of K^CrgAsg, we obtain a set of param¬ 
eters (in eV) as following: £4 = 0.2274, £4 = —0.0573, 
t 2 = -0.0040, t\ = -0.0481, £4 = 0.0242, t' 2 = 0.0057, 
£3 = 0.0147, £5 = -0.0023, £3 = 0.0022, £ 5 = -0.0124, 


£' 3 = 0.0015, £5 = -0.0058, £3 = -0.0020, t' 5 = -0.0297. 
We also obtain the chemical potential p c hem = 0.2365 
eV. As shown in Fig. [21 the set of parameters well repro¬ 
duces the Fermi surfaces. They also reproduce the energy 
dispersion along high symmetry lines and the density of 
states (DOS) near the Fermi level well 111. 


FIG. 2: Comparison of the Fermi surfaces for K 2 CrzAsz ob¬ 
tained from the DFT calculation (a), and from the present 
model Hamiltonian (b). The a, /3 and 7 bands are labelled. 
The T point is at the center of the Brillouin zone. 



By diagonalizing Hq we obtain three energy bands with 
dispersion and eigen-wavefunctions <f>k„, hereafter p 
is the band index. The bare susceptibility tensor xo(q) 
is defined as 


mn.m n / \ 

x 0 (q) 


-to 

k /iv 



m' * 
k/i 


^k+q// 


.ra* 

k+qzy 


^k fi ^k+qiy 


x[/(e kAt ) - /(e k+q*/)], (5) 


where f(e k^) is the Fermion occupation number, N is the 
number of unit cells. Note that Xo(q) can be also written 
as a 9 x 9 matrix in the two-body orbital space spanned 
by the basis (|m) ( 8 > |?r) : m, n = 1 , 2 ,3}. 


Interaction: Now we consider the electron-electron in¬ 
teraction. In the spirit of Hubbard approximation, i.e. 
retaining only the intra-unit-cell terms, we obtain the in¬ 
teracting Hamiltonian which respects D^h symmetry, 11] 
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Hint — 2 “I 'y^, El IV 1 C Lj C Lj C ‘™ C ‘™ + '^ c imo- C lm5- C *"to' C *”icr + J {c\ma C \ma c i3crCi3a + /l.C.)] 

2 f 771—1 £7 
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+ El y^(^2Cj TOa .c| ) ^ cr ,CiTOcr'Ci TOO - + •^c| ma .cj ) ^ a .,Ci mcr 'CiTOcr + t^2 C imcr C i3CT' Ci 3cr' c ima- + J c i m cr C i3a-' C ™<r' C i3cr) 


m= 1 crcr' 


+ El ^l c ]’3o- C i3g C »3°' C *3o-1* , 


( 6 ) 


where m is the opposite orbital to m within E' repre¬ 
sentation (1 = 2 and 2 = 1), and er is the opposite 
spin to er. Here U\ , U 2 and J are intra-orbital repul¬ 
sion, inter-orbital repulsion and Hund’s coupling for the 
two E' states, which satisfy that U\ = U 2 +2 J. U[ is the 
intra-orbital repulsion for the A\ state, and J' involves 
one of E' states and the A\ state. 

By Fourier transformation, we have 

= E E r r’ mV (k,k-q;k',k' + q) 

kk'q mnm'n' aa' 

Xc kma C k'm / a ,Ck '+ c l n ' a '' Ck ~cincn (7) 

where the bare vertex function r™ n ’ m n (k, k — q; k', k' + 
q) can be written as a matrix in the two-body orbital 
space, 

f 0 (k, k - q; k', k' + q) = (1 - <W)f 8 + <Wf t , ( 8 ) 


where f* is bare vetex in the equal-spin channel and P s is 
in the opposite-spin channel. The two-body orbital space 
can be decomposited to subspaces {| 11), |22), |33)} © 
{|12), |21}} © {113), 113)} © (|23), |32)}. Thus, f s and 
ft are block diagonal in this set of basis, 

/ U x U 2 U-2 \ 

f s = u 2 u x u’ 2 © jn© j'n © j'n, (9a) 

\V 2 u: 2 U[J 

( o u 2 ui\ 

ft = \ U 2 0 U' 2 © Jcri © J'fTi © JVi, (9b) 

V [/' f/' 0 J 


where rri is Pauli matrix and n = 1 + <ti isa2x2 matrix. 

Effective pairing interaction. For weak coupling, the 
full vertex function ” (k, k — q; k', k' + q) can be 

evaluated diagrammatically, for instance, through the 
random phase approximation (RPA). To study the SC 
pairing, we only need to keep the vertices in pairing chan¬ 
nels with k' = —k, say, f CTCT /(k, k'; —k, —k') = 14,7 (k, k'). 
Therefore, tE(k,k') = (1 - <W)H s (k, k ') + <WVt(k, k') 
serves as an effective pairing interaction to study super¬ 
conductivity instability. 

In study of the cuprates, Scalapino et al. used RPA to 
calculate V aiy (k,k') in a single Hubbard model. 12] The 


RPA involves two types of Feynman diagrams in addition 
to the bare vertex function Vo(k, k') 11]. One contains 
the bubble diagrams, and the other contains the ladder 
diagrams with Cooperon. Here we generalize the calcu¬ 
lation in Ref. [3 to the multi-orbital case. The effective 
pairing interaction I4(t)(k, k') from the bubble diagrams 
is 


E 6 $(k,k') = i{(f t + f 8 )[i + xo(q)(ft+f 8 )] 


-1 


T(f t - f 8 )[l + xo(q)(f t - f,)]” 1 ]^) 


where © takes — for s and + for t. The effective pairing 
from the ladder diagrams is 


fs(t)(k> k 7 ) = f s(t) xo(p)f s(t) [l - Mp)fs(t)}-\ (11) 

where p = k + k' and q = k — k' [ll| . For the notation, 
matrix A is related to matrix A via the following relation, 

A , , — A , , 


We can also project effective pairing potential V aCT / (k, k') 
into three single particle bands through 


fT'(k,k')= £ 


±m* rm * 


y k/, V'— 


mn,m n 


(k,k')C^- k v ; 


mnm' n' 

. (12) 

where cf>^ is the single particle eigenwavefunction in the 
/ 7 -band. 

Superconducting pairing instability. Because the three 
Fermi surfaces have different shapes, k and — k are always 
in the same band for k near the Fermi surfaces. For weak 
coupling, we only consider intra-band pairing (k to — k). 
A single band gap function can be written as 


A (k) = * [croV’(k) + cr • d (k)] cr 2 , (13) 

where (k) = if) (—k) is the spin-singlet gap function, 
and the d-vector d (k) = — d (—k) describes the spin- 
triplet pairing, ao is the unit matrix, er 1 , 2,3 are Pauli ma- 
trice. To measure the intrji-band SC pairing instability, 
we follow Scalapino et cd. .fl2(| and introduce a dimension¬ 
less coupling constant 

/ St f S3 (km (k, k ')<c (ko 
" W/SE^I^k )! 2 
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where f dS is the integration over the p -band Fermi 
surface and is the Fermi velocity, the /x-band form 
factor ^ CT ,(k) oc A cr( 7 /(k). 

In the absence of the spin-orbit coupling, the spin- 
singlet component 1 / (k) will not mix with the spin-triplet 
component d (k) in Eo. lfTdl) . The possible single-band SC 
gap functions on hexagonal lattice are listed in Table U 
up to the first and second NN bonds [l3|. 

TABLE I: Superconducting gap functions -0(k) and d(k) on 
hexagonal lattice. The gap functions are intra-band and are 
classified according to D^h group irreducible representation 
T. The functions s x 2 +y 2 , p x , p y , d x 2 _ y 2 , and d xy have been 
defined in Eq. ©■ Other functions are defined as follows, p z = 
sin k c , d z 2 = cos k c , f x ( x 2 - 3 V 2 ) = sin k a +sin kb — sin(fc a + kb), 

fy(3x 2 -y 2 ) = sin (2fc a + k b ) - sin (fc a + 2k b ) - sin ( ka - k b ). 


r 

spin-singlet V/k) 

spin-triplet d(k) 

a; 

1? $x 2 -\-y 2 5 d z 2 

fy(3x 2 -y 2 )Z 

a 2 


fx(x 2 —3 y 2 )Z 

E' 

(d x 2 — y 2 , d X y) 

(p x ,p v )z,p z (x,y) 

A'/ 

Pz fy(3x 2 —y 2 ) 

PxX + p y y,p z z 

A'/ 

Pz fx(x 2 — 3y 2 ) 

PyX - p x y 

E" 

Pz(Px,Py) 

(p x x-p y y,p y x + p x y) 


We next examine all the pairing states in Table Q] to in¬ 
vestigate which pairing channel will dominate. For sim¬ 
plicity, we set U\ = U[ = U, U' 2 = Ui and J' = J at 
first. The relevant A^’s are plotted as functions of J/U 
in Fig. [3] and those A M not plotted have either negative 
or negligibly small values. 

Fig. 3(a) shows the pair coupling constant in vari¬ 
ous channels for U = 0.5 eV. Results for U < 0.5eV 
are similar. The spin-triplet f y ( 3 X 2 - y 2 ) state in the 7 - 
band appears at a finite J/U and become dominant when 
J/U > 1/3. The f y ( 3 x 2 -y 2 ) state has a line nodal gap. 
Fig. 3(c) plots the results at 17 = 2 eV, to represent large 
U. In that case, the spin-triplet p z z state at the a-band 
dominates in all the realistic region of J/U < 0.45. The 
p z z state has a full gap at the quasi-ID a-band Fermi 
surface. In the intermediate region, e.g., U = l.OeV, 
p z z and fy( 3 x 2 - y 2 ) state compete against each other, p z z 
dominates at small J/U, while f y ( 3 X 2 - y 2 ) state become 
strong at large J/U. We have also calculated the pair 
coupling constant for J' ^ J while keep U\ = U[ = U 
and U 2 = U 2 - The results are similar in a wide param¬ 
eter region 0.5 < J'/J < 2.0. Note that at large U and 
J/U, the SC pairing is found in the /3-band with pairing 
symmetry p x x ± p v y and p y x ± p x y , which gives point 
nodal in the gap function. However, for these point 
nodally gapped states are tiny. 

Note that the divation of our tight-binding model from 
DFT results in some details would not change the above 
statements qualitatively. Because it is the Fermi surface 
shape and DOS near the Fermi level that determine 
and pairing symmetry in weak coupling. 



J/U 

FIG. 3: The dimensionless pair coupling constant A M vs 
J/U for (a) U = 0.5eV, (b) U = l.OeV and (c) U = 2.0eV. 
Here we have set Ui = U[ = U , U 2 = U 2 and J' = J. All 
the dominated states are spin-triplet pairing states, (a) U = 
0.5eV, f y ($ x 2 _ y 2 ) appear and dominate when J/U > 1/3. (b) 
U = l.OeV, p z z state and f y ( 3x 2 _ y 2 ) compete against each 
other, (c) U = 2 . 0 eV, p z z state dominates over most J/U. 


Discussions and Summary We would like to mention 
some issues which have not been addressed in this paper 
but may be interesting for future work, (i) We neglect 
spin-orbit coupling in the theory for simplicity. Actually 
the D^h lattice breaks inversion symmetry, and the spin- 
orbit coupling may mix spin singlet and triplet states 
within the same D^h irreducible representation in Table 
m (ii) We have only considered intra-unit-cell interaction 
in the Hubbard model. Since the Wannier functions are 
extended in a CrgAsg cluster, the inter-unit-cell interac¬ 
tion along the c-axis could be sizable, (iii) Since molec¬ 
ular orbitals are more extended than atomic orbitals, U 
and J are estimated to be smaller than or comparable to 
the bandwidth. In this paper we start with weak cou¬ 
pling to study SC instability. However, A 2 Cr 3 As 3 lattice 
is quasi-ID, an alternative approach would be to model 
the system as a coupled Tomonaga-Luttinger liquid. 

In summary, we have proposed a minimal model 
to study superconductivity in A 2 Cr 3 As 3 (A=K,Rb,Cs), 
which involves three molecular orbital states in each unit 
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cell. With the help of symmetry, we have deduced a 
tight-binding model with 3 molecular orbitals for the sys¬ 
tem, which compares well with the results of the density 
functional theory. We have derived effective pairing in¬ 
teractions within the RPA, and found that the dominant 
pairing channel is always spin-triplet. For small U, a 
spin-triplet state with line nodal gap, / y ( 3 a , 2 _ y 2 ), at the 
3D 7 -band will dominate at moderate Hund’s coupling. 
While for large U, a spin-triplet fully gapped state, p z z, 
will dominate at the quasi-ID a-band. The state we 
find at small U appears to be most relevant to the com¬ 
pound. The pairing state of f y ^ x I 2 - y 2 ) has nodal lines on 
planes k a = kb, k a = —2kb and kb = —2k a , say, T-K-L-A 
planes on hexagonal lattices. The T-K-L-A planes cross 
small sections at y-Fermi surface. Our theory at small 
U predicts line zeroes in gap function and appears to be 
consistent with existing experiments showing non-BCS 
gap function and particularly the low temperature Lon¬ 
don penetration depth measurement. Our prediction can 
be tested in further experiments including angle resolved 
photoemission spectroscopy. 
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Supplementary Material: Theory for superconductivity in alkali chromium arsenides A2Cr 3 As 3 

(A=K,Rb,Cs) 


D 3h GROUP: CHARACTER TABLE, CLEBSCH-GORDAN COEFFICIENTS AND OTHERS 


In this section, we provide useful informations on the point group D 3 h, which have been used to determine the 
formatoin of Hamiltonian and analyze superconducting pairing symmetry. 


Character table 


D 3 h group is a direct product of D 3 h group and a Zi group generated by the horizontal reflection ah, namely, 
L> 3 h = -D 3 g) ah ■ The character table for D 3 h is the following. 

Below we list some functions for D 3 h irreducible representations on hexagonal lattice, limited to the second NN 
bonds. 
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TABLE II: Character table for D 3 h group. 


D 3 h = D 3 x ah{6m2 ) 

E 

2 C 3 

3C 2 

<7h 

2 S3 

3(7 v 

x 2 + y 2 , z 2 

x{x 2 - 3 y 2 ) 

A'\ 

1 

1 

1 

1 

1 

1 


y(3x 2 - y 2 ) 

a ' 2 

1 

1 

-1 

1 

1 

-1 

(x 2 ~y 2 ,xy) 

(• x,y),(xz 2 ,yz 2 ) 

E' 

2 

-1 

0 

2 

-1 

0 

yz(3x 2 - y 2 ) 


A" 

1 

1 

1 

-1 

-1 

-1 

xz(x 2 — 3 y 2 ) 

z,z 3 

AS 

1 

1 

-1 

-1 

-1 

1 

(xz, yz) 

(z(x 2 - y 2 ),zxy) 

E" 

2 

-1 

0 

-2 

1 

0 


TABLE III: Functions for D 3 h group representation on hexagonal lattice. 


D 3h = D 3 x a h (6m2) 

COS (ka + kb) + COS ka + COS kb 

sin k a + sin kb — sin (k a + kb) 

a:[ 


sin (2 k a + k b ) — sin (k a + 2 kb) — sin (fc 0 — kb) 

a' 2 

(cos k a + cos kb — 2 cos (k a + kb), \/3(cos kb — cos k a )) 

(2 sin(fc a + kb) + sin k a + sin kb, -\/3(sin kb — sin k a )) 

E' 



A" 


sin k c 

a' 2 ' 

sin k c { cos k a + cos kb — 2 cos (k a + kb), %/3(cos kb — cos k a )) 

sin k c {2 sin(fc a + kb) + sin k a + sin kb, %/3(sin kb — sin k a )) 

E" 


Since 


when feo <C 1 , we have 


k 


a 


h 


-k x - —k 
2 2 

-k x + —k 
2 2 


y 


cos k a + cos kb — 2 cos (k a + kb) 

•\/3(cos kb — cos k a ) 
2 sin(/c a + kb) + sin k a + sin kb 
-\/3(sin kb — sin k a ) 


\{k 2 x - k 2 y ), 

--k k 
3k x , 

Sky. 


Clebsch-Gordan coefficients 

We shall also find out Clebsch-Gordan (CG) coefficients for later use. Since D 3 h = D 3 ® ah, we only need to 
consider the CG coefficents for D 3 group. For D 3 group, there are three representations A\. A 2 and E, where Ai and 
A 2 are ID and E is 2D. The natural notations for A\ and A 2 are 0 and 0 respectively, say, |0) and | 6 ) denote Ai and 
A 2 states respectively. We also denote the two E states as |1) and | — 1). Here |1) and | —1) are two eigenstates whose 
eigenvalue under C 3 rotation are and e~ 1 ^ respectively. The two E states can be also written as |x) and | y), 
which are related to |± 1 ) as 

\±l) = ^=(\x)±i\y)). 

For the decomposition Ti ® T 2 = ®F, we have 

|0) = |0,0) = 10,0> = —(|1,—1) + | —1,1)), 

|o) = |o,o) = |o,o> = -^= (|i,—1) - |-i,i)), 

|i) = |o, 1} = 1 1,0) = |o, 1) = 1 1, o> = |—1, —1), 
l-i) = |o,-i> = |-i,o) = - |o,-i) = - |-i,o) = |i,i>. 
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In terms of |x) and | y), we also have 


| 0 ) 

|o> 

|x) 

I y) 


1 

7! 

i 

V2i 


(\x,x) + | y,y )), 

: i\x,y) - | y,x )), 


|0, x) = \x, 0) = |6,x) = \x, 5) = -j={\x,x) - | y,y)), 


-1 


10 , y) = \y, 0 } = 16, y) = \y,0) = -^=(| x,y) + \y,x}). 


The CG coefficients for group can be extended to D^h group straightforward, via the following correspondence 
between D^h irreducible representations and £>3 irreducible representations, 


A" — (ThA'i Al, 

A '2 = <7hA' 2 A 2 , 
E" = a h E'^E. 


How x, y, z transfer under ah 


To analyze the spin-triplet pairing states, one need to study how the basis x. y. z transfer under D^h symmetry 
operations, especially under the operation cry- The basis states x,y,z are defined as 

x = |sr = o> = 7=Htt> + IW>), 
v = \ S T = o) = 7= (ltt> +III)), 

£ = \sr = 0 > = 7=(in} + i^t)). 

Note that the horizontal reflection can be written as 

Z7T 

CTh = IC 2z = I exp(-—er 2 ) = -Ha z , 

where / is the inversion operation. Since angular momentum operators and their eigenstates will not change under I, 
we have 


Oh It) = 
oh II) = * II) • 


Therefore 


a h x 


OhV 


OhZ 


Itt) + III)) = —x, 
~^2 Jh d)')') HD) = — & 
—f^Oh (III) + lit)) = Z. 


We would like to thank Daniel F. Agterberg for pointing out this correct transformation, who found a similar table 
as Table I in the main text. 


DERIVATION OF Ho AND COMPARISON WITH DFT CALCULATION 

A generic three-orbital tightbinding model without spin-orbit coupling can be written as 

Ho = c kms£kmnCkns- 

kmns 



In the basis {|1}, |2), |3)}, the Hermitian matrix £k can be expressed in terms of Gell-Mann matrices, 


Ik = ^2 ££A T 


An — 


A 3 — 


An — 


\ 

i cO 


/ 0 -i 0 \ 

,Ai = 

1 0 0 

) A 2 = 

i 0 0 

1 

v o 0 0 j 


O 

O 

O 


T —0 

where GellMann matrices Ai_s are defined with the unit matrix Ao as follows, 

/ 1 0 0 
0 1 0 
v o 0 1 

/ 1 0 0 
0-10 
v 0 0 0 

/ 0 0 0 
0 0 1 
v o 1 0 



1 0 0 1 N 


^0 0 -i\ 

,A 4 = 

0 0 0 

5 A 5 = 

0 0 0 

/ 

l 1 0 V 


v i 0 0 J 


, A7 = 


0 


( 0 0 
0 0 

v 0 i 0 


, — \ — 



The issue is how to determine Writing ^ = (cj^ , c]^, cj^ ), then we can find out how ¥>k s A r </?ks t rans f er 
under D^h symmetry operations. From the CG coefficients, we find that they can be classified to different irreducible 
representations of D$h group. 

A^ • ^ks^O^ksi ^ks'^^ks 

A 2 : ^’k s A 2 </?ks; 

E' ■■ (^LAiV^ks, ¥>L A 3^k s ), (<^LA 4 ^k s , ^L A 6^k s ), (¥>LA5¥> ks , ^LArV’ks). 


Since Hq itself is invariant under any D^h symmetry operation, thereby belongs to A[ representations, we require 
that ££ belong to different irreducible representations as follows, 

A : && 

A' 2 : & 

E’ : &,&),(&&),(£,&)■ 

Considering A' 2 representation need longer bonds and (£k,£k) w iU break time reversal symmetry in general, we have 
£k = ^k = sk = 0- With the help of CG coefficients, we obtain Eq.(4) in the main text. 

Fig ISlI shows that the set of tight-binding paramters given in the main text well reproduce the energy dispersion 
along high symmetry lines and the density of states near the Fermi level. 


DERIVE THE THREE-ORBITAL HUBBARD MODEL 


The electron field operator ipa(f) can be expanded in terms of a complete set of Wannier functions, 

ipa{r) = EE Wim (X)Cim<j i 

m i 

where Ci ma annihilates an electron with orbital m and spin a at lattice site i. A generic interacting Hamiltonian is 
given by 

Hint = ^E / Xa(n)X a '{r 2 )V(r 1 -r 2 )'ipa'{r2)'4’a{n)df 1 df2 : 


where V ex t(r) is external periodic potential, and V{f\ — r 2 ) is the screened Coulomb interaction. Hq can be written as 
a tight-binding model and we shall focus on the interacting part H. int . In the spirit of a Hubbard type approximation, 
i.e. retaining only the terms on the same lattice site, we have 


^ 1 Vmm',n'n c \ m o C i m ' (T > c in , <T' c in0') 


h ^ = \Y. e 

i mnm'n' aa 


j j 
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K-points 

(a)Energy dispersion 



FIG. SI: Comparison of the energy dispersion along high symmetry lines (a) and the density of states (b) for K 2 CrsAs 3 
obtained from the DFT calculation and from the present model Hamiltonian. The red and blue lines are the results obtained 
with DFT calculation and the model calculation, respectively. 


where 


Vmm',n'n = J (f 2 )V (fi - f 2 )w in ' {f 2 )w in (fi )df 1 df 2 . 

Let’s consider nonvanishing V mm ' :n ' n according to D^h group symmetry. Since the interaction V{f\—f 2 ) = V(|fi — f 2 \) 
respects the full point group symmetry and belong to the A\ representation. Wether a V mm ',n'n vanishes can be 
determined by the CG coefficients. 

Firstly, for to, n = 1, 2, the nonzero can be written in terms of the integrals explicitly, 


mm,mm 


= J \u)im(ri)\ 2 V(r! - f 2 )\w im (f 2 )\ 2 df 1 df 2 = U 1 , 

Vmm,mm = J \ Wim (bl )\ 2 V {fi - f 2 )\vJifh (f 2 ) | 2 df=i df 2 = U 2 , 

Vmm,mm = Jw* m (fi)Wim(ri)V(fi-f 2 )w* m (f 2 )w irn (f 2 )dfidr 2 = J, 
= J w* m (fi ) Wim{ri)V(ri - r 2 )w m {r 2 )w* m (f 2 )dr\dr 2 = J*. 


mm,mm 


where U\ is the intra-orbital interaction and U 2 is the inter-orbital interaction. If one chooses the Wannier function 
to be real, then J* = J, and J is the Hund’s exchange energy. 

Secondly, we can choose that |1) transfers as x and |2) transfers as y under D^h symmetry operations. Note that 
V{fi — r 2 ) is invariant not only under D^h symmetry operations, but also under all the 0(3) symmetry operations. 
So that V xx ,xx = Vyy,yy = U\ ■ Under a O 3 rotation along the c-axis, the two Wannier functions | 1 ) and \2) (denoted 
by x and y) transfer as 


w x cos 6w x + sin 9w y , 
w y —> — sin 9iu x + cos 9w y , 

where 9 = -£■. The integral V XXtXX should keep invariant under this operation. Assuming w x and w y are real, after 
straightforward algebra, we have 

U\ = (cos 4 9 + sin 4 9)U\ + 2U 2 cos 2 9 sin 2 9 + 4 cos 2 9 sin 2 9 J, 

where the transfermation peroperties under C 2 are used. From the above, we have 


U\ = U 2 + 2 J. 




















10 


Thirdly, we involve the state |3). For to, to' = 1,2 (or x, y), the relevant nonvanishing terms are given in the 
following, 

V 33,33 = J\w i3 (ri)\ 2 V (n - r 2 ) | Wi 3 (r *2 )| 2 dridr 2 = U[, 

V 3 mi m3 = J \wirn{r 1 )\ 2 V(r 1 -r 2 )\wi 3 {r 2 )\ 2 dridr -2 = U' 2l 

V m3 ,m3 = j w* m (ri)wi3(fi)V (fi - r 2 )w* 3 {r2)wi m (r2)dr 1 dr2 = J ', 

Vmm',33 = J Wi m (n)wi 3 (ri)V(f! - f 2 )w^ ml (r 2 )wi 3 (f 2 )df 1 df 2 = J" m -. 

When Wannier functions are real, = J'■ For m ^ to', we set to = x and m' = y, thus under the C 3 rotation, 


J xx = cos 2 0J" X + sin 2 dJyy + 2 cos 6 sin 9J xy = J xx + 2 cos 0 sin 9J xy . 


The above gives rise to 


J"y = 0. 


In summary, we have 


Hint — „ ^ El (Ul c \mcr C lmcr' C i™<7' C iTn<r + ^2 c lrocr C lmcr' c im<r'Cjma + J C imcr'Cjmcr J c \ma C \ r na' C imiT'Cjma) 


i.aa' \ m—1 


+ ^l c i3o- C l3cr ,Ci 3o- ,c i3o- + El[^2 c Imo- C ls3o-'^sBcr'Gmcr + J cJ ma .cj 3a .,Ci mC r'Ci3cr + J i,c\ miJ c\ ma ,Ci3a'Ci3 a + Il.C.)] > 

m J 

x 2 

= JEE W^ima^imS^maCima + J c \m<T C ima C imv c im<j + J ( c jma C ima C i3a c i3a + ^-C.)] 


2 ,cr m=l 

l 


A 9 El El (^ C lmaCifha' Cirncr 1 Cirna + J Eicr^roer' 0 *”^' Cimcr + V2 c i rna - c i3a' c i3a' Cima + J ^ima^sa' C inaa’ c i3a) 


i.cra' m =1 


; El ^l C l3cr C i3CT C *3a-Cj3<7; 


where to is the opposite orbital to to within the E representation, and a is the opposite spin to a. 


RPA FOR EFFECTIVE PAIRING INTERACTION 


We shall generalized the calculation by Scalapino et al. to the multi-orbital case. Diagramatically, RPA contains 
two types Feynman diagrams beside the bare vertex function Vo(k, k') as shown in Fig. [32] One contains the bubble 
diagrams, the other contains the ladder diagrams with Cooperon. 

For the single orbital case, the RPA effective pairing interaction for spin-singlet and spin-triplet channels are 


Vs(k, k') = 


U 


_ U 2 x oCk + kQ 

l-t/ 2 y§(k-k') l-t/xo(k + k')' 

U 2 xo(k-k') 


Ft(k,k,) 

Note that Fig. IS2T c) is absent in U s and Fig. ISfil' bf. (d), (e) and (f) are absent in V t because a' = d in the U term. 
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FIG. S 2 : Feynman diagram for RPA calculation of effective pairing interaction V (k, k'). (a) Effective vertex function V (k, k'). 
(b) Bare vertex function Vo(k, k'). (c) and (d) Bubble diagrams, (e) and (f) Ladder digarams with Cooperon. 


For multi-obtial case, we shall replace the numbers U and xo by the bare vertex functions r s ( t ) in Eq. (9) and 
susceptibility tensor xo in Eq. (5) in the main text. The effective pairing interaction E s ( t )(k, k') from the bubble 
diagrams are 

V s bub (k,k') = f s - f s xo(q)f t - f t xo(q)f s 

+fsXo(q)ftxo(q)f* + f*xo(q)f s xo(q)ft + ftxo(q)ftxo(q)fs + f s xo(q)f s xo(q)f s + • • • 
+r s xo(p)r s + f s xo(p)r s xo(p)r s d— 

= 2 + r s )[l + Xo(q)(f t + f s )] -1 - -(f t - f s )[l + X' 0 (q)(ft - fs)] -1 


and 


vr 6 (k,k') = r t - r s xo(q)r s - r t xo(q)r t 

+ftXo(q)ft%o(q)ft + ftxo(q)f s xo(q)fs + f s xo(q)fao(q)fs + fsXo(q)fsXo(q)ft + • 
+ftXo(p)ft +ftXo(p)ftXo(p)ft H- 


— -(ft + f s )[i + x°(q)(ft + fs)] 1 + — f s) [i + xo(q)(ft — f. 

while those from ladder diagrams are 

Vf d (k, k') = f sXo (p)f s + fsXo(p)fsXo(p)fs + • 


s;j > 


= fsXo(p)fs[l - Xo(p)fs]" 


and 


v t lad (k,k') = r t xo(p)r t + r t xo(p)r t xo(p)r t + • • ■ 

= f*Xo(p)ft[l — Xo(p)ft] -1 ) 

where p = k + k' and q = k — k', and the matrix A is related to matrix A through the following relation, 

A , , — A , , 

mn.m'n' — ^-nn' ,m'm- 
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Matrix summation in the bubble diagrams 

Below we write down the details on summing over the two series of matrice in the bubble diagrams, 

A = f s - f s xo(q)f t - f t xo(q)f s 

+r s xo(q)r t xo(q)r t + f t xo(q)f s xo(q)f 4 + f t x 0 (q)ftXo(q)fs + f s x 0 (q)f s x 0 (q)f s h— 


and 


B = f t - f s xo(q)f s - f t xo(q)r t 

+ftXo(q)ftXo(q)f* + ftXo(q)f s xo(q)f s + f«xo(q)ftxo(q)f« + f s xo(q)f«xo(q)ft h— 

where every term in A cotains odd number of F s and every term in B contains even number of F s , and each xo 
contribute one minus sign. To do this summation, we multiple Xo(q) at both sides at first, then add them together 
or substract one to another. Flence 

xo(q)s + xo(q)i = xo(q)(f t + f s )[i+ xo(q)(f t + f s )] _1 , 
xo(q)£ - xo(q)A = xo(q)(f t - f s )[i+ xo(q)(f t - f s )] _1 . 

Then we obtain that 

A = i(f t +f 5 )[l+xo(q)(ft+f s )]- 1 -i(ft-f s )[l + xo(q)(ft-f s )]- 1 , 

B = i(f t +f s )[l+xo(q)(f t +f s r 1 + i(f t -f s )[l + xo(q)(f t -f s )]- 1 . 

INTEGRATION OVER THE FERMI SURFACE 

To calculate the dimensionless coupling constant A, we need to integrate over the Fermi surface. However, it is 
difficult to do this kind integration numerically. Insteadly, we recall the method to calculation DOS. The DOS is 
derived as follows, 

f dS'k <5fc(k) _ r dS'k 1 
[e, ~J s{e) ( 2tt)3 de - J (27t) 3 |v k |" 

Numerically, we use the following Lorentz lineshape to calculate it, 

D{ - e) = Ev 2 (e k — e) 2 + 5 2 ’ 

k 

where N is the number of fc-points in the first Brillouin zone and 6 is the width of Lorentz lineshape. Thus, the 
integration in A can be converted to 

^Ek.k' ^J^+si ^MViKk.^gCk') 

^E k (6k _ M V+^ lg( k )l 2 



